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Abstract 

The covariant density functional theory with the point-couphng interaction PC-PKl is compared 
with new and accurate experimental masses in the element range from 50 to 91. The experimen- 
tal data are from a mass measurement performed with the storage ring mass spectrometry at 
Gesellschaft fiir Schwerionenforschung (GSI) [Chen et al, Nucl. Phys. A 882 71 (2012)] . Al- 
though the microscopic theory contains only 11 parameters, it agrees well with the experimental 
data. The comparison is characterized by a rms deviation of 0.859 MeV. For even-even nuclei, 
the theory agrees within about 600 keV. Larger deviations are observed in this comparison for the 
odd-A and odd-odd nuclei. Improvements and possible reasons for the deviations are discussed in 
this contribution as well. 

PACS numbers: 21.10.Dr, 21.60.Jz, 21.30.Fe 
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Nuclear masses play a vital role not only in the nuclear physics but also in weak- 
interaction studies and astrophysics. In particular, the mass of heavy neutron-rich nuclei 
is one of the basic topics in nuclear physics and is essential for understanding the rapid 
neutron-capture nucleosynthesis process (r process). 

During recent decades, great achievements in mass measurements of neutron-rich nuclei 
have been made, thanks to the applications of cyclotron, storage ring, and penning trap 
facilities Despite the experimental progress, the masses of a large number of neutron- 
rich nuclei relevant to the r process remain unmeasured due to difficulties in production, 
separation, and detection. Therefore, reliable theoretical predictions of the nuclear mass are 
imperative at the present time. 

The theoretical determination of nuclear masses can be traced to the von Weizsacker 
mass formula, which was proposed based on the famous liquid drop model (LDM) ^. 
Since then, tremendous effort has been made in pursuing different possible extensions of 
the LDM which are known as the macroscopic-microscopic models, for example, the finite- 
range droplet model (FRDM) [s], the extended Thomas- Fermi plus Strutinsky integral with 
shell quenching (ETFSI-Q) model J], and the Weizsacker- Skyrme (WS) model 5|. It turns 
out that these macroscopic-microscopic models work pretty well in the description of known 
nuclides, but their predictions show a large deviation for very neutron-rich nuclides. In 



addition, the local mass relations such as the isobaric multiplet mass equation (IMME) 



8| are 



the Garvey-Kelson (GK) relations and the residual proton-neutron interactions 
also used to give predictions of unmeasured masses. 

In principle, an ideal mass formula would be one in which the masses of all nuclei are 
derived from the basic nucleonic interactions. In this regard, the microscopic-rooted mass 
model, which treats the macroscopic part and the microscopic corrections in a unified frame- 
work, is usually believed to have a more reliable extrapolation to the unknown regions. In 
the past decade, a series of microscopic-rooted mass models based on the nonrelativistic den- 
sity functional theory (DFT) have been developed with the Hartree-Fock-Bogoliubov (HFB) 
method and have achieved great success in the description of known masses (see Refs. 9l-lll| 
and references therein). In these models, the model parameters are fitted to essentially all 
the experimental mass data. 

Apart from the nonrelativistic DFT, the covariant density functional theory (CDFT) 
has also received wide attention due to its successful description of many nuclear phenom- 
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15|. There exist a number of attractive features in the CDFT, especially in its 
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practical applications in the self-consistent relativistic mean-field (RMF) framework. The 
most obvious one is the natural inclusion of the spin degree of freedom, and the relativistic 
effects are responsible for the existence of the approximate pseudospin symmetry 
in the nuclear single-particle spectra and spin symmetry in the antinucleon spectra 
Moreover, it is of particular importance that the CDFT includes nuclear magnetism 
which play s an important role in the microscopic description of the nuclear magnetic mo 
ments 
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28| and nuclear rotations [29|-|32|. 



The first RMF mass table was reported in Ref. [33j for about 2000 even-even nuclei 
with 8 < Z < 120 up to the proton and neutron drip lines but without including pairing 
correlations. Later on, the ground-state properties of 1315 even-even nuclei with 10 < Z < 



98 were calculated by including the pairing correlations with BCS method [3J|. In Ref. 35| . 
by using a state-dependent BCS method with zero-range 6 force, the first systematic study 
of the ground-state properties of more than 7000 nuclei ranging from the proton drip line 
to the neutron drip line was performed with the meson-exchange effective interaction TMA. 



This mass table works well in the r-process nucleosynthesis calculations 36 



Very recently, a new point-coupling effective interaction PC-PKl has been proposed by 
fitting to observables of 60 selecte d sp herical nuclei, including the binding energies, charge 



radii, and empirical pairing gaps 



40|. This effective interaction particularly improves the 



description for isospin dependence of binding energies and it has been successfully used in 
describing the Coulomb displacement energies between mirror nuclei 4]J], fission barriers 42|, 



nuclear rotations 
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J-|32|, etc. On the other hand, a new opportunity has been opened for a 
crucial test of the predictive power of the theories with the comparison to accurate new mass 
values of 53 heavy neutron-rich isotopes from Sn to Pa measured with the storage ring mass 
spectrometry at GSI 43|. This mass measurement is characterized by a small systematic 
error of about 10 keV, which is valid for all data. The overall mean experimental accuracy 
is about 19 keV. From these 53 values, 31 have been measured for the first time, whereas 
for the additional 22 nuclides the experimental error has been significantly improved. Since 
none of these data have been used to determine the parameters of the CDFT with PC-PKl 
and most of them have not been used in the fitting of the widely used mass models in the 
market, these data provide a very good test for the mass-prediction power for nuclei in this 
region. The work for all the measured masses with CDFT definitely should be done in the 
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future. 

In this work, the CDFT with the point-couphng interaction PC-PKl is apphed to investi- 
gate the new experimental masses of 53 heavy neutron-rich isotopes at the storage ring mass 
spectrometry at GSI, including 31 cases measured for the first time 43|. The theoretical 
predictions are compared with the experimental data as well as other theoretical results. 

The starting point of the CDFT is a general effective La gra ngian density where the 



nucleons are coupled with zero-range point-coupling interaction |40l. |44|. |45| . By means of the 
mean-field approximation and the no-sea approximation, one could obtain the corresponding 
Kohn-Sham equation, which has the form of a Dirac equation 

[a-p + (3{m + S) + V]iJk = ek^k, (1) 
with the scalar S{r) and vector V{r) potentials, 

S{r) = asps + (3spI + IspI + SsAps, (2a) 
V{r) = avpv + IvPv + <^yApy + eA^ 

+aTvT3PTV + Stv^TsPtv- (2b) 

Here, m is the nucleon mass, and as, ay, dTV, Ps^ Is, Tv? ^S: ^v, and Stv are the coupling 
constants. The Coulomb field Aq is determined by Poisson's equation. 

The iterative solution of these equations yields the single-particle energies and expectation 
values of total energy, quadrupole moments, etc. Here, we present only the expectation value 
of total energy 



= j d^r V-^i-(^) (" ■ P + Mr) 

, "5 2 , 3 , 75 4 , "^S A 

+ + yPs + -jPs + yPsAps 
+yPy + yPy + yPvApy 

+^^Ptv + ^PTV^PTV + ^eAop^v | > (3) 

where vl is the occupation probabilities of single-particle state and psir), pv{f), and Ptv(^) 
are the local densities in scalar, vector, and isovector-vector channels, respectively. 

For open-shell nuclei, pairing correlations are taken into account by the BCS method 
with a zero-range 6 force. Thus, we have to add to the functional Eq. ([3]) a pairing energy 



term, 

E,.. = -J2 T I d'r^*rir)KAr), (4) 



4 

T=n,p 



where Vr is the constant pairing strength, and the pairing tensor ^(r) reads 

{r) = -2j2fkUkVk\Mr)\^ (5) 



k>0 



with /fc being a smooth cutoff factor |46|, |4j 

As the translational symmetry is broken in the mean-field approximation, the center-of- 



mass (cm.) correction shon 
microscopic cm. correction 



d be made for the spurious cm. motion. Here, we adopt the 
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with A being the mass number and Pc.m. = '^f Pi being the total momentum in the cm. 
frame. 

Similarly, the rotational symmetry is also violated for the deformed nuclei. Therefore, we 
further introduce the correction energy 

i?rot = ~^{J') (7) 
with J the angular momentum operator and X the moment of inertia calculated from the 



Inglis-Belyaev formula 



50|-|52| 



Finally, the total energy for the nuclear system reads 

E = Edf + E'pair + Ec.m. + -E'rot, (8) 

and the corresponding binding energy is Eb = —E. 

In this work, the Dirac equation is solved on the basis of an axially deformed harmonic 



oscillator potential [53|. A basis of 16 major oscillator shells is used in the calculations, and 
convergence has been tested in the calculations with 18 major shells. 

In Table [H the calculated binding energies E^'^^, the rotational correction energies (RCEs) 
= —Ej-ot, and the quadrupole deformations P2 for nuclei with masses measured at the 



storage ring mass spectrometry at GSI 43[ are listed together with the differences between 



the data E^^^ and the calculated binding energies AE = E^^^ — E^^^. 

It is found that most nuclei are deformed except the nuclei with proton number close to 
the magic numbers of 50 and 82. Note that it is not necessary to make rotational correction 



to the binding energy for a spherical nucleus, since it preserves the rotational symmetry. 
Therefore, we here consider the RCEs only for the deformed nuclei with \j32\ > 0.02. The 
present calculations with PC-PKl reproduce the experimental data very well and the root 
mean square (rms) deviation results in 0.859 MeV for 53 nuclei, and 0.805 MeV for the 31 
cases measured for the first time. 

TABLE 1: The experimental binding energies E^^^, 
calculated binding energies E^^^, rotational correc- 
tion energies E^^, and quadrupole deformations (32 
with the PC-PKl, as well as the differences between 
the data and the calculated results AE. The root 
mean square (rms) deviation A is listed in the last 
row. The energies are shown in million electron volts 
(MeV). 



Element 


Z 


A 






^B 


132 


AE 


Sn 


50 


128 


1077.20 


1078.67 


0.00 


0.00 


-1.47 


Sb 


51 


133 


1112.28 


1113.70 


0.00 


-0.02 


-1.42 


Te 


52 


136 


1131.28 


1131.48 


1.15 


0.09 


-0.20 


La 


57 


144 


1192.29 


1191.36 


0.71 


0.19 


0.93 


Ce 


58 


146 


1208.38 


1208.25 


2.00 


0.20 


0.13 


Pt 


78 


202 


1591.54 


1590.09 


0.00 


0.00 


1.45 


Au 


79 


202 


1592.40 


1590.78 


0.57 


-0.08 


1.62 


Hg 


80 


207 


1624.09 


1624.17 


0.00 


-0.02 


-0.08 


Tl 


81 


213 


1653.45 


1652.30 


0.00 


0.00 


1.15 


Bi 


83 


217 


1677.18 


1677.10 


0.87 


-0.04 


0.08 


Bi 


83 


218 


1680.77 


1679.28 


0.04 


0.07 


1.49 


Po 


84 


219 


1688.58 


1687.34 


0.33 


0.11 


1.24 


Po 


84 


220 


1694.07 


1694.04 


1.44 


0.10 


0.03 


Po 


84 


221 


1697.62 


1696.21 


0.34 


0.12 


1.41 


Po 


84 


222 


1702.98 


1702.97 


1.54 


0.11 


0.01 
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TABLE I: (Continued). 



Element 


Z 


A 


rnExp 




^B 


/32 


AE 


At 


85 


220 


1694.15 


1694.33 


0.45 


0.12 


-0.18 


At 


85 


221 


1699.82 


1701.19 


1.47 


0.12 


-1.37 


At 


85 


222 


1703.72 


1703.84 


0.55 


0.13 


-0.12 


At 


85 


223 


1709.31 


1710.26 


1.21 


0.13 


-0.95 


At 


85 


224 


1713.10 


1712.53 


0.17 


0.14 


0.57 


Rn 


86 


223 


1711.56 


1711.44 


0.67 


0.15 


0.12 


Rn 


86 


224 


1717.56 


1718.35 


1.59 


0.15 


-0.79 


Rn 


86 


225 


1721.56 


1720.80 


0.41 


0.16 


0.76 


Rn 


86 


226 


1727.40 


1728.04 


1.71 


0.16 


-0.64 


Rn 


86 


227 


1731.32 


1730.65 


0.81 


0.18 


0.67 


Rn 


86 


228 


1737.06 


1737.41 


1.80 


0.17 


-0.35 


Ft 


87 


224 


1717.41 


1716.91 


0.42 


0.14 


0.50 


Fr 


87 


225 


1723.46 


1724.75 


1.79 


0.14 


-1.29 


Fr 


87 


226 


1727.81 


1727.27 


0.64 


0.15 


0.54 


Fr 


87 


227 


1733.74 


1734.87 


1.90 


0.16 


-1.13 


Fr 


87 


228 


1738.13 


1737.28 


0.64 


0.17 


0.85 


Fr 


87 


229 


1743.89 


1744.20 


1.93 


0.20 


-0.31 


Fr 


87 


230 


1748.12 


1747.22 


0.92 


0.18 


0.90 


Fr 


87 


231 


1753.64 


1753.79 


2.05 


0.21 


-0.15 


Ra 


88 


231 


1756.69 


1756.10 


1.81 


0.24 


0.59 


Ra 


88 


232 


1762.48 


1762.41 


2.27 


0.24 


0.07 


Ra 


88 


233 


1766.72 


1765.72 


1.69 


0.26 


1.00 


Ra 


88 


234 


1772.23 


1772.06 


2.26 


0.25 


0.16 


Ac 


89 


229 


1747.27 


1747.36 


1.88 


0.23 


-0.09 


Ac 


89 


230 


1752.19 


1751.21 


1.12 


0.24 


0.98 


Ac 


89 


231 


1758.34 


1758.28 


1.90 


0.24 


0.06 


Ac 


89 


232 


1763.02 


1762.50 


0.66 


0.26 


0.52 
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TABLE I: (Continued). 



Element 


Z 


A 








/32 




Ac 


89 


233 


1768.94 


1768.69 


1.90 


0.25 


0.25 


Ac 


89 


234 


1773.47 


1772.41 


1.52 


0.27 


1.06 


Ac 


89 


235 


1779.03 


1778.66 


1.88 


0.26 


0.37 


Ac 


89 


236 


1783.24 


1782.32 


1.52 


0.27 


0.92 


Th 


90 


235 


1781.57 


1780.76 


1.80 


0.28 


0.81 


Th 


90 


236 


1787.40 


1787.27 


2.27 


0.27 


0.13 


Th 


90 


237 


1791.77 


1791.08 


1.84 


0.28 


0.69 


Pa 


91 


235 


1782.49 


1782.02 


1.40 


0.27 


0.47 


Pa 


91 


236 


1787.52 


1786.25 


1.12 


0.28 


1.27 


Pa 


91 


237 


1793.40 


1792.37 


1.03 


0.27 


1.03 


Pa 


91 


238 


1798.10 


1796.41 


0.73 


0.28 


1.69 



A 0.859 



In order to present the calculated results more plainly, in Fig. [H we show the differences 
between the experimental data 43| and the calculated binding energies with the PC-PKl [4^ 
in comparison with those obtained from the CDFT calculations with TMA 35| as well as 
the mass models HFB-21 [uj and FRDM [31 . 

It is found that the deviations given by both PC-PKl and TMA are less than 1 MeV for 
most nuclei. Both cases underestimate the binding energies for nuclei near N = 126 and 
overestimate them for nuclei near = 82. Moreover, the deviations for nuclei ^^^La and 
^■^^Ce in TMA are even as large as 3 MeV. This indicates that PC-PKl is more robust in the 
predictions of nuclear masses, especially for the neutron-rich nuclei. The up-to-date nonrela- 



tivistic mass model HFB-21 



ll| could also reproduce the data quite well. However, it should 



be noted that there are more than 20 parameters in HFB-21 including the phenomenologi- 
cal corrections, and they are determined by constraining the nuclear-matter parameters and 
optimizing the fit to the full data set of the 2149 measured masses with N, Z > 8 (both 
spherical and deformed) [5J|, including some of the present data (improved masses in the 
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work 43|), while the functional PC-PKl contains only 11 parameters, which are determined 
by fitting to the binding energies of only 60 selected spherical nuclei and charge radii for 17 
ones 40|. It should be mentioned that the macroscopic- microscopic FRDM also gives good 



agreement with the data. 
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FIG. 1: (Color online) Differences between the experimental data [43] and the calculated binding 
energies by the covariant density functional theory (CDFT) with the effective interaction PC- 



PKl 



40| and TMA 



35l | as well as those from the mass models HFB-21 



111] and FRDM 



a. 



Focusing on the heavy mass region in Fig. [H especially for the nuclei with A > 210, one 
can clearly see some systematic regularities for the mass deviations. Both the PC-PKl and 
TMA results show linear-like increasing tendency with the mass number, and the HFB-21 
results exhibit the opposite tendency. In addition, the results given by PC-PKl could be 
roughly classified into two branches. For FRDM, all the deviations for the heavy mass region 
are concentrated around 0.6 MeV. 

In order to investigate this in more detail, in Fig. [21 differences between the data and the 
calculated binding energies for nine isotopic chains from Z = 83 to Z = 91 are shown. It 
can be seen that the deviations for FRDM remain constant with mass number. For both 
PC-PKl and TMA, the binding energy differences of most isotopic chains grow steadily 
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with the neutron number, which indicates that the isotopes in most isotopic chains become 
more underbound with the increase of the neutron number. In contrast, for the mass model 
HFB-21, the isotopes in most isotopic chains become more overbound with the increase of 
the neutron number. Such difference may result from the different recipes employed in the 
calculations for treating the pairing correlations. In the CDFT calculations with PC-PKl 
and TMA, the pairing correlation is treated with the simple BCS method, while in the 
mass model HFB-21, it is treated with the sophisticated Bogoliubov transformation. For 
the neutron-rich nuclei, the latter could describe the pairing correlation with the continuum 
properly and therefore include more correlations which makes the nucleus more bound. 
It was found in Ref. 55| that there is a general tendency for overbinding with the HFB 
calculations in comparison with the HF -|- BCS case, and the corresponding mass deviations 
could be up to 1.5 MeV for open-shell nuclei. Therefore, future efforts should be devoted 



to the sophisticated and heavy Bo goliu bov calculations with continuum for t 



rich heavy nuclei as shown in Refs. 56|, l57l| for the spherical case and Refs. 58 
deformed case. 
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FIG. 2: (Color online) Same as Fig. [T]but for the isotopic chains from Z = 83 to Z = 91. 
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Comparing the results given by PC-PKl and TMA in Fig. [2l one could see that there 
are noticeable odd-even staggering in the results given by PC-PKl, which eventually leads 
to the two linear-like branches of the deviations. This staggering phenomenon influences 
the rms value to a great extent. In Fig. [31 the rms deviations of the calculated binding 
energies by the CDFT with PC-PKl and TMA as well as those from the mass models HFB- 
21 and FRDM are shown for the even-even, odd-A, odd-odd, and the whole set of nuclei, 
respectively. For the whole set of nuclei, the rms deviation given by PC-PKl is about 0.8 
MeV. This is larger than the rms value 0.7 MeV given by HFB-21 but smaller than the value 
of 1.0 MeV given by TMA. Speciflcally, for the 12 even-even nuclei, PC-PKl achieves very 
good agreement with the data and the rms deviation is about 0.6 MeV. This is at the same 
good level as the HFB-21 mass model and much better than that given by TMA. However, 
for the 25 odd-^ nuclei and 16 odd-odd nuclei, the rms values given by PC-PKl become 
larger but still within 1 MeV. In contrast to the microscopic frameworks, the mass model 
FRDM gives very good agreement for the odd-odd nuclei with a rms deviation less than 0.5 
MeV. This lowers the corresponding rms deviation for the whole set of nuclei to 0.6 MeV. 




Even-Even Odd-A Odd-Odd Total 

FIG. 3: (Color online) The rms deviations of the calculated binding energies by the CDFT with 
the effective interaction PC-PKl and TMA as well as those from the mass models HFB-21 and 
FRDM for the even-even, odd-A, odd-odd, and the whole set of nuclei. 

To explore the reason for this enhancement of the rms values for odd particle systems, in 
Fig. m correlations between the deviations of the calculated binding energies with PC-PKl 
from the data and the calculated rotational correction energies are shown for even-even, odd- 
A, and odd-odd nuclei respectively. It is clear that the rotational correction energy plays 

11 



an important role in the determination of the nuclear masses with PC-PKl. In particular, 
for the odd-odd nuclei, the rotation correction energies are systemically smaller than those 
of the even-even nuclei even if they have similar deformations. Apparently, this leads to the 
systematical underestimation of the binding energies for odd-odd nuclei. 




Even-Even Odd-A Odd-Odd 

FIG. 4: (Color online) Correlations between the deviations of the calculated binding energies with 
PC-PKl from the data and the calculated rotational correction energies for even-even, odd-A, and 
odd-odd nuclei. Note that the order in different data sets is independent. 



This may result from two reasons. On one hand, for the odd particle system, the odd 
nucleon breaks the time-reversal symmetry and the time-odd components of the vector 
potential do not vanish. Although it is generally believed that the contributions due to 
the time-reversal symmetry breaking contribute little to the bulk properties such as binding 
energies, they may significantly influence features like magnetic moment 2J] or moment 
of inertia 6ll|, which depend on a few valence nucleons. Here, the calculation of rotation 
correction energy with Eq. (JT]) involves with the moment of inertia, which is calculated 
without the time-odd components. Therefore, the neglect of the time-odd components may 
affect the rotation correction energy. 

On the other hand, for odd particle systems, the unpaired particle will block its occu- 
pied level in the BCS calculations; that is the Pauli principle prevents this level from the 
scattering process of nucleon pairs by the pairing correlations. In principle, the ground- 
state for odd particle systems should be the state with the lowest energy determined by 
calculations blocking each possible level near the Fermi surface. However, this is incredibly 
time-consuming. 

Here the blocking method adopted in the BCS method is the following: For each step of 
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the self-consistent calculation, the level to be blocked is determined by the filling approx- 
imation from the single-particle spectra obtained. Of course, it may occur that the state 
chosen by this procedure is not the true ground state. As pointed out in Ref. |62|, such 
wrong selection of the blocking state generally leads to a difference less than 0.2 MeV for 
the binding energy. However, this would remarkably influence the moment of inertia and 
thus the rotation correction energy. Therefore, further improvement in the evaluation of 
the rotation correction energy, in particular the moments of inertia, would be needed in the 
framework of CDFT. 

In summary, the CDFT with the point-coupling interaction PC-PKl is applied to investi- 
gate the new experimental masses of 53 heavy neutron-rich isotopes in the element range of 
50 to 91. The experimental data are from a mass measurement performed with the storage 
ring mass spectrometry at GSI 43|]. The functional PC-PKl contains only 11 parameters. 



which are determined by fltting to the binding energies of 60 spherical nuclei and charge 
radii for 17 ones. It is found that the CDFT with PC-PKl can reproduce the experimental 
data quite well and the corresponding rms deviation is 0.859 MeV. For the 25 odd-A nuclei 
and 16 odd-odd nuclei, the rms values given by PC-PKl are a little large but still within 1 
MeV, which may be improved in the future by treating properly the time-odd component 
and the moment of inertia. An excellent predictive power with an accuracy of about 600 
keV has been achieved for even-even nuclei with the PC-PKl. 
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